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ALsirnct 

Wingtip  Pnd  fuselage  mounted  side  force  generating  (SFG)  surfaces 
were  designed  end  installed  on  the  XBQM-106  remotely  piloted  vehicle  (RPV) 
to  enhance  its  lateral  terminal  response  characteristics.  These  surfaces 
uere  sized  and  positioned  in  an  attempt  to  keep  the  net  rolling  and  yawing 
moments  about  the  CG  unchanged  when  the  aircraft  vas  side  slipping.  The 
7LEXSTAB  digital  computer  system  in  conjunction  with  traditional  hand 
calculated  methods  ’-'ere  used  to  evaluate  the  REV's  stability,  control,  and 
time  response  characteristics.  Because  of  the  propeller  location  and  the 
unconventional  shape  of  the  fuselage,  three  mathematical  models  were  gener- 
ated for  comparing  lateral  characteristics:  a.  basic  fuselage  and  vertical 
tail  effectiveness  unchanged,  b.  forward  fuselage  effectiveness  reduced 
by  50%.  c.  vertical  tail  effectiveness  increased  by  50%.  Comparisons  of 
the  three  data  sets  with,  the  hand  calculated  results  indicated  that  the  model 
with  increased  vertical  tail  effectiveness  more  closely  modeled  the  unmodified 
airplane.  In  addition  the  simulation  showed  that  the  initial  sizing  and 
locations  of  the  SFG  r>roduced  rolling  and  yawing  moments  under  side  slip 
conditions.  The  time  history  responses  to  rudder,  aileron,  rudder/aileron, 
and  wind  gust  inputs  ware  generated  and  plotted  by  the  FLEXSTAB  program. 

From  these  plots  it  v.’C  determined  that  the  SFG  modification  increased  the 
yaw  rate  response  10.6'»  while  decreasing  the  roll  rate  10.7^  and  the  side 
.:lip  angle  response  25.''%.  The  dutch  roll  damping  was  increased  A%  end  the 
DR  period  decreased  The  roll  mode  time  to  half  amplitude  increased  16.7^ 

and  the  spiral  stability  increased  for  the  modified  vehicle.  The  longitudinal 
stability  suffered  a slight  decrease  due  to  the  addition  of  the  SFG. 
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DESIGN  AND  EVALUATION  OF  A SIDE  FORCE  GENERATOR 
( 

MODIFICATION  FOR  THE  XBQM-106,  REMOTELY  PILOTED  VEHICLE 


I . Introduction 

Background 

After  WWII  our  military  strategies  began  an  increasing  shift  toward 
nuclear  warfare  and  as  a result  the  role  of  conventional  weaponry  was 
downgraded.  The  recent  Vietnam  war  and  Middle  Eastern  conflicts  have 
shown  that  conventional  warfare  is  still  a very  important  aspect  of 
general  warfare.  In  light  of  this  fact  and  our  increasing  technology 
the  Air  Force  determined  in  1971  that  there  exists  a need  for  an 
effective,  low  cost,  expendable  weapon  system  capable  of  being  remotely 
or  automatically  deployed  for  use  against  both  stationary  and  mobile 
ground  targets.  The  Air  Force  Flight  Dynamics  Laboratory  (AFFDL)  at 
Wright-Patterson  AFB,  OH  was  given  the  task  of  providing  a Remotely 
Piloted  Vehicle  (RPV)  to  serve  as  a test  bed  vehicle  to  evaluate  the 
electronic  seekers  and  sensors  that  will  be  utilized  by  this  weapon 
system.  This  study  deals  with  the  latest  of  these  RPV's,  the  XBQM-106 
(Fig  1). 

Statement  of  the  Problem 

The  problem  can  be  defined  in  this  manner:  To  design  a side  force 
generator  modification  for  the  XBQM-106  and  evaluate  its  effect  on  the 
vehicle's  response  and  stability.  Side  force  generators  were  studied 
because  this  modification  should  improve  the  lateral  response  character- 
istics, in  sped  Fie,  the  yawing  response  capability. 
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Figure  1.  XBQM-106  Basic  Model 


Objectives 

The  primary  obj'ective's  of  this  report  are: 

1.  Prepare  design  drawings  of  suitable  side  force  generating  (SFG) 
surfaces  and  submit  them  to  AFFDL  for  construction  and  instal- 
lation on  the  XBQM-106. 

2.  Calculate  the  stability  and  control  derivatives  for  both  the 
basic  and  modified  vehicles. 

3.  Calculate  the  dynamic  response  of  both  vehicles  to  given  inputs 

4.  Determine  the  change  in  stability  characteristics  resulting 
from  the  SFG  modification. 

Scope 

This  study  describes  the  stepwise  process  of  sizing  the  SFG  system 
for  the  XBOM-106.  Stability  and  control  derivatives  for  both  vehicles 


are  calculated  using  the  digital  computer  program  FLEXSTAB  (Refs  1,2,3) 
and  partially  validated  by  using  traditional  geometrically  based  param- 


eter estimation  techniques.  With  this  data  the  basic  and  modified 
models  were  compared  using  computer  generated  time  history  responses. 


Assumptions 

FLEXSTAB,  the  primary  analytical  tool  used  in  this  thesis,  is 
based  on  linear  aerodynamic  theory  and  consequently  the  user  must  avoid 
such  conditions  that  would  produce  flow  seperation  and  its  resulting 
nonlinearities . 

Due  to  the  complexity  of  the  flow  associated  with  the  propeller  and 
other  similar  problems  with  modeling  the  propeller,  a torque  free  envir- 
onment excluding  power  effects  and  "corkscrew"  propeller  wash  was  assumed 
in  this  study. 

All  aeroelastic  effects  are  assumed  to  be  negligible  and  therefore 
all  the  analyses  in  this  report  are  based  on  the  rigid  airplane. 

In  hand  calculating  the  derivatives  it  was  asstimed  that  interference 
effects  between  different  aircraft  bodies  (for  example,  horizontal  tail/ 
vertical  tail)  were  negligible.  The  FLEXSTAB  programs  do  calculate 
interference  effects  so  this  assumption  does  no.t  apply  in  the  computer 
studies . 

The  assumed  coordinate  system  centered  at  the  center  of  gravity  (CG) 
is  depicted  in  figure  2. 
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Figure  2.  Coordinate  System 
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Approach 

To  begin  this  evaluation  the  side  force  generators  were  designed 
for  the  lower  Kinj.tips  and  fuselage  bottom.  The  upper  SFG  was  then 
sized  to  counter  balance  the  rolling  moment  due  to  yaw  caused  by  the 
first  two  surfaces. 

After  the  gv^cMetry  of  the  modified  vehicle  was  established  the 
problem  of  para;nelcr  estimation  was  confronted.  Hand  calculated  esti- 
mation methods  f.'ou  several  sources  were  used  to  determine  some  of  the 
more  important  s::;ibility  and  control  derivatives. 

Next  both  vehicles  were  defined  mathematically  aind  evaluated  by 
the  FLEXSTAB  cor.ipuLcr  program.  The  previously  mentioned  hand  calculated 
derivatives  were  used  as  a compari:-.on  with  the  computer  generated 
numbers . 

FLEXSTAB  was  then  used  to  generate  static  and  dynamic  stability 
characteristics  for  each  model.  From  these  analyses  time  histories 
were  produced  using  various  control  and  wind  gust  inputs  as  forcing 
functions.  The  conclusions  conceri’.ing  the  effect  of  the  SFG  modification 
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II.  Design  of  Side  Force  Generators 


Introduction 

This  section  contains  a detailed  explanation  of  the  placement  and 
design/sizing  process  used  to  locate  the  side  force  generators  (SFG)  on 
the  vehicle.  Symmetric  airfoils  (NACA  0015)  were  chosen  for  use  as  the 
SFG  because  of  strength  reasons  and  good  performance  in  the  low  subsonic 
mach  number  range,  .04  to  .3. 

Placement 

The  first  step  in  the  design  process  was  the  placement  of  the  SFG. 
The  goal  was  to  locate  the  surfaces  in  such  a manner  that  horizontally 
(x  direction)  the  centers  of  pressure  would  coincide  with  the  fore/aft 
center  of  gravity  location  and  vertically  (z  direction)  so  that  the  SFG 


2.  In  conjunction  with  number  1,  the  end  plate  effect  decreases 

the  induced  drag  and  increases  the  lift-curve  slope,  . 

a 

These  benefits  are  particularly  important  for  performance 
considerations. 

3.  Referencing  Figure  3,  one  can  see  the  ailerons  lie  adjacent  to 
the  wingtips.  In  this  region  the  pressure  distribution  is 
improved  and  consequently  the  aileron  control  power,  C , is 
increased. 

Because  the  CG  is  located  below  the  mainwing  (see  Figure  6)  a 
balancing  SFG  had  to  be  placed  on  the  lower  section  of  the  fuselage. 


Main  Wing 


Figure  6.  Center  of  Gravity  Location 


Figure  7.  Aerodynamic  Side 
Forces  on  SFG  caused  by 
Sideslip 
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Placing  the  SFG  here  provides  a force  below  the  CG  which  counterbalances 
the  forces  acting  above  the  CG  (see  Figure  7).  Forces  A and  B create  a 
moment  in  a counterclockwise  direction  and  C provides  a balancing 
moment  in  the  clockwise  direction. 

Sizing 

The  actual  sizing  process  was  constrained  in  two  ways.  The  maxi- 
mum vertical  span  under  the  wingtip  was  restricted  to  7.3".  This 
allows  the  RPV  a maximum  of  ±10*  of  bank  when  the  center  skid  is  in 
contact  with  the  ground.  It  was  felt  that  this  amount  of  roll  capability 
would  be. needed  to  assist  in  lateral  dir'ctional  control  during  landing 
phases.  The  other  constraint  was  that  tj:e  lower  fuselage  SFG  had  to  be 
stowed  in  the  fuselage  daring  the  launcl. i .'jg  sequence.  This  required 
that  the  span  of  this  surface  be  apprcxi  .ately  12".  Considering  these 
restrictions,  the  planforms  of  the  lower  .vingtip  SFG  and  fuselage  SFG 
were  designed  as  shown  in  Figures  8 and  " , respectively. 
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with  the  placement  and  sizing  of  tiiese  two  SFGs  complete  the 
rolling  moment  due  to  side  slip  produced  by  them  had  to  be  calculated. 
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Figure  9.  Fuselage  SFG 


A seven  step  process  was  used  to  evaluate  the  rolling  moment  component 
of  each  SFG. 

1.  Compute  planform  area,  S: 

S = 3^  b (C^  + . (1) 

where  b is  the  span,  C is  the  tip  chord  and  C is  the  root  chord. 

T R 

Note:  All  of  the  side  force  generators  were  trapezoidal  shapes. 

2.  Calculate  aspect  ratio,  JR: 

JR  = /S  (2) 


3.  Calculate  the  lift-curve  slope,  C^^  : 


1 + a / /S 
o 


(Ref  4:4) 


' (3) 


where  a is  the  theoretical  section  lift-curve  slope  (assumed  to  be  2-n 
o 

in  these  calculations).  Since  side  force  was  needed  rather  than  lift, 

it  was  assumed  that  C % C , where  C is  the  coefficient  of  side 

yg  'a 

force  due  to  a change  in  side  slip  angle. 


4.  Calculate  taper  ratio,  X: 

X = C^/Cj^  (4) 


5.  Determine  the  loeation  of  the  center  of  pressure,  CP: 


l-h2X 

3(1+X) 


(Ref  5:502) 


(5) 


6.  Calculate  the  resulting  side  force,  F , caused  by  side  slip 


angle: 


F 

- y 


3 q s 


(6) 


where  q is  the  dynamic  pressure. 

7.  Determine  the  rolling  moment,  L,  about  the  CG  created  by  the 
force  acting  through  its  CP. 

L = F X t (7) 

y 

where  I is  the  distance  the  CP  is  above  or  below  the  CG. 

Knowing  the  net  rolling  moment  caused  by  the  lower  wingtip  SFG 
and  fuselage  SFG,  the  dimensions  of  the  upper  SFG  were  changed  in  an 
iterative  manner  to  create  an  equal,  opposing  moment.  The  planform 
of  this  surface  is  depicted  in  Figure  10. 
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H;>nd  Calculated  Parzimcter  Estimation 


Introduction 

This  section  lists  in  tabular  form  various  stability  and  control 
derivatives  that  were  hand  calculated  and  the  reference  used  to  obtain 
the  method.  The  express  purpose  for  these  calculations  was  to  obtain 
approximate  values  to  be  used  in  validating  the  FLEXSTAB  model.  Also 
included  in  this  section  is  the  method  by  which  the  center  of  gravity 
and  moments/crossproduct  of  inertia  were  calculated. 


Hand  Calculated  Derivatives 

Since  there  exists  no  wind  tunnel  or  flight  test  data  for  the 
XBQM-106,  the  basic  geometry  (see  Table  I)  was  the  only  input  used  for 
determining  the  derivatives. 


Table  I 

Geometry  Information 


Mainwing 

Horizontal  Tail 

Vertical  Tail 

Planform  Area  (in) 

2714 

523 

470.4 

Span  (in) 

139 

48 

24 

Mac  c (in) 

19.52 

11 

19.6 

Aspect  Ratio 

7.12 

4.364 

1.224 

Taper  Ratio 

.751 

.692 

.352 

Sweep  Angle  - LE 

2.5* 

4.5“ 

44'“ 

Sweep  Angle  - C/4 

1.0® 

2.5“ 

39.5“ 

Incidence  - Root 

4* 

0“ 

- 

Incidence  - Tip 

1“ 

0“ 

Twist 

3“ 

0“ 

- 

Dihedral 

0“ 

0“ 

Several  sources  were  consulted  during  the  search  for  parameter  estima- 
tion methods.  The  (wo  primary  sources  used  in  this  study  were: 


I 

I 

I 


i 


"Stability  and  Control  of  Airplanes  and  Helicopters,"  (Ref  4),  and 
"Methods  for  Estimating.  Stability  and  Control  Derivatives  of  Conven- 
tional Subsonic  Airplanes,"  (Ref  6). 

A standard  cruise  flight  condition  within  the  vehicle's  operating 
envelope  was  arbitrarily  chosen  and  used  in  the  calculations,  see 
Table  II. 

Table  II 


Flight  Condition 


Altitude  . 

2000'  MSL 

Forward  Velocity 

110.9  ft/sec 

Mach  Number 

0.1 

Weight 

140  lb  Basic/145  lb  SFG 

Flight  Path  Angle 

0°  (Straight  S Level) 

Bank  Angle 

0“  - 

All  the  calculations  were  based  on  a power-off  situation  conse- 
quently no  account  was  taken  for  such  unknowns  as  gyroscopic  force  of 
the  propeller,  torque  of  the  engine,  and  "corkscrew"  prop  wash  that 
the  tail  section  experienced..  This  assumed  condition  was  not  intended 
to  infer  that  the  above  mentioned  items  are  negligible,  but  that  their 
incorporation  was  beyond  the  scope  of  this  study. 

Listed  in  Tables  III  and  IV  are  the  derivatives  and  the  references 
from  which  the  methods  of  calculation  were  taken, 

A certain  amount  of  caution  needs  to  be  exercised  when  using  these 
numbers  because  it  is  extremely  difficult  to  arrive  at  accurate 
solutions  using  only  the  geometry  of  the  vehicle.  The  following  quote 
was  taken  from  Ref  6. 
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Table  III 


Hand  Calculated  Longitudinal  Derivatives 


Derivative 

Value  (rad  ^) 

Reference (s) 

S 

o 

.0162* 

4:165 

c. 

4.933 

4:5,  6:3.2, 

L 

a 

- 

7:460 

S 

a 

.2939 

6:3.1 

-1.115 

6:3.3,8:3.40 

6.4916 

6:5.1 

s 

-17. .1028 

6:5.2 

■ 

.7448 

6:10.3 

C 

-2.5153 

6:10.4 

* Dimensionless 


"When  better  accuracy  is  desired,  it  is  always 
recommended  to  use  (wind)  tunnel  data  and/or  a 
combination  of  (wind)  tunnel  data  and  the 
methods  of  Reference  1 (Datcom)." 


In  addition  to  the  uncertainties  associated  with  a particular 
method,  the  values  obtained  for  a given  derivative  using  different 
methods  can  also  vary  significantly.  As  stated  previously  the  purpose 
of  generating  these  hand  calculated  derivatives  was  only  to  provide 
a basic  check  on  the  FLEXSTAB  model. 
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Table  IV 

Hand  Calculated  Lateral  Derivatives 


Derivative 

Value  (rad 

Refercnce(s) 

C 

-.4333 

4:223,  6:7.1 

-.0621 

6:7.3 

3 

C 

.068 

4:239 

c 

-.011 

6:8.1 

p 

-.452 

4:231 

c 

n 

P 

.0052 

6:8.3 

C 

.3713 

6:9.2 

.1608 

4:222,  6:9.3 

C 

n 

-.2177 

6:9.3 

r 

^9. 

6A 

.2458 

6:11.1 

C 

-.0252 

6:11.2 

C 

y^R 

.1971 

6:12.1 

6R 

.0017 

6:12.1 

C 

"6R 

-.1203 

6:12.1 

Inertia/Center  of  Gravity  Calculation 


The  moments/cross  product  of  inertia  and  the  CG  were  calculated 


using  the  digital  computer.  The  process  consisted  of  dividing  the 


vehicle  into  several  finite  volumetric  elements,  calculating  the  weight 


of  each  element,  and  calculating  its  distance  from  an  arbitrary 


reference  point.  Using  this  information  the  program  solved  for 

I , I , I and  the  center  of  gravity  location.  Sec  Table  V. 
yy  7.7.  xz  H j 


I 


I 


■f 

4' 

1 


I 


t'l 

{' 


S 

1. 
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.Table  V 

Moments/Crossproduct  of  Inertia 


Basic  Model  (sl-ft  ) 

SKG  Model  (sl-ft  ) 

^xx 

9.167 

13.125 

^yy 

25.486 

25.625 

32.431 

36.25 

^xz 

-1.528 

-1.59 

Weight 

140  lb 

145  lb 

Note  the  large  increase  (43%)  in  the  second  moment  about  the  x 
axis,  Ijjjj.  See  Appendix  A for  a listing  of  Inertia/CG  Program. 
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IV.  FLEXSTAB  Parameter  Estimation 

Introduction 

This  section  explains  the  applicable  FLEXSTAB  programs  and  how  they 
were  used  to  predict  the  stability  and  control  derivatives  and  model 
information  of  the  XBQM-106. 

FLEXSTAB  Explanation 

The  FLEXSTAB  2.01.00  Computer  Program  system  (FLEXSTAB)  is  a system 
of  digital  computer  programs  based  on  linear  aerodynamic  theories  for 
evaluating: 

(1)  Static  stability 

(2)  Dynamic  stability 

(3)  Trim  states 

(4)  Linear  systems  analysis 

(5)  Structural  loading 

(6)  Elastic  deformation 

The  vehicle  being  analyzed  can  be  of  an  arbitrary  shape/configuration 
in  subsonic  or  supersonic  flow. 

The  system  consists  of  fourteen  programs,  eight  of  which  deal  with 
airplane  definition  and  the  remaining  six  are  airplane  analysis  programs. 
In  this  study  only  five  of  the  programs  were  utilized.  They  are  listed 
in  the  order  used:  Geometry  Definition  program  (GD),  Aerodynamic 
Influence  Coefficient  program  (AIC),  AIC  Matrix  Correction  program  (CAIC), 
Stability  Derivatives  and  Static  Stability  program  (SDSSS)  and  Time 
Histories  program  (TH).  See  Appendix  B for  the  general  FLEXSTAB  overview. 


1 . 

I! 

ii 

I . 


if 
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Geometry  Overivew 

The  sole  input  to  this  program  is  the  vehicle's  dimensions  which 

were  obtained  from  the  original  design  blueprints.  The  program  takes 

this  input  and  creates  a three-dimensional  mathematical  geometry 
definition  of  the  airplane.  The  basic  model  is  shown  in  Figure  11  and 
the  model  with  the  SFG  modification  is  depicted  in  Figure  12.  The 

reader  should  note  the  increased  panel  density  in  the  proximity  of  the 

wing  tips  of  Figure  12.  This  was  done  to  increase  the  accuracy  of  the 
modeling  process  in  these  regions  of  changing  pressure.  This  geometry 
definition  consists  of: 

(1)  Slender  bodies  - A mathematical  model  used  to  define  the 
fuselage  and  nacelles.  Slender  bodies  are  bodies  of  revolution 

with  an  area  distribution  similar  to  the  area  distribution  of 
the  part  it  is  representing. 

(2)  Thin  bodies  - A mathematical  model  used  to  define  wing,  hori- 
zontal and  vertical  tails.  Thin  bodies  are  defined  by  a 
planform  that  is  identical  to  the  actual  airplane's  planform. 

(3)  Interference  bodies  - A mathematical  model  used  to  model  the 
interference  effects  between  a slender  body  and  its  neighboring 
thin  bodies.  It  is  composed  of  a shell  that  surrounds  a 
slender  body  in  regions  where  there  are  significant  wing/body 
Interference  effects# 

Due  to  the  type  of  analysis  in  this  study,  thin  bodies  were  given 
no  camber  or  thickness.  This  is  a reasonable  approach  since  incremental 
effects  of  the  modification  were  the  main  interests  in  the  study. 

Due  to  its  unconventional  shape,  the  forward  fuselage  had  to  be 
modeled  as  a thin  body  connecting  two  slender  bodies.  The  problems 
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Figure  12.  SFG  FLEXSTAB  Model 


2 


associated  with  this  type  of  modeling  will  be  discussed  later.  See 
Appendix  B for  a flow  chart  of  the  Geometry  Definition  program. 

Aerodynamic  Influence  Coefficient  Program 

This  portion  of  FLEXSTAB  uses  the  mathematical  geometry  definition 
(GD)  and  computes  the  steady  and  unsteady  matrices  used  to  define  the 
airplane  aerodynamically . Basically  the  AIC  program  uses  a modification 
of  linear  aerodynamic  theory  introduced  by  Woodward  (Ref  9)  to  determine 
the  effect  of  each  panel /control  point  in  the  geometry  on  every  other 
panel/control  point.  This  is  an  extremely  long  program  in  terms  of 
computer  time.  For  example,  typical  durations  are:  3000  seconds  compila- 
tion time  and  2500  seconds  input/output  time.  See  Appendix  B for  a flow 
chart  of  the  AIC  program. 

AIC  Matrix  Correction 

CAIC  is  a program  that  allows  the  user  to  modify  an  existing  AIC 
matrix.  In  this  study  the  forward  vertical  portion  of  the  fuselage 
was  modeled  using  a thin  body.  The  data  gathered  using  this  scheme 
suggested  that  this  thin  body  was  more  effective  than  the  actual  fuselage 
mainly  because  the  magnitude  of  the  C term  was  very  small.  For  this 

reason  the  CAIC  program  was  utilized  to  reduce  the  effect  of  this  Thin 

body  by  50%.  A second  modification  of  the  AIC  matrix  was  accomplished 

by  increasing  the  effect  of  the  vertical  tail  by  50%.  To  differentiate 

between  the  data  gathered  using  the  three  different  AIC  matrices,  the 

following  convention  was  established.  ".5C  " denotes  data  where  the 

^FUS 

front  fuselage  effect  was  reduced  by  50%.  "Baseline"  represents  data 

from  the  unaltered  AIC  matrix  and  "1.5 G " denotes  data  from  a model 

, ^VT 
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with  a 50%  increase  in  the  effect  of  the  vertical  tail.  A more  detailed 


explanation  of  the  reason  for  modifying  the  AIC  matrix  is  presented  in 
the  SD  C SS  section.  See  Appendix  B for  a flow  chart  of  the  CAIC 
program. 


Stability  Derivatives  and  Static  Stability  Program 

This  program  uses  the  geometry  definition  and  the  AIC  matrix  to 
compute  rigid  static  and  dynamic  stability  derivatives  and  mode  shape 
characteristics  of  the  model. 

Tables  VI  and  VII  list  the  longitudinal  and  lateral  derivatives 
respectively  for  both  the  basic  and  modified  vehicles  using  the  baseline 
AIC  matrices.  Also  for  comparison  purposes,  the  hand  calculated  values 
of  selected  derivatives  are  included  in  these  tables.  Note;  These 
values  are  based  on  the  flight  condition  in  Table  II. 

Comparing  the  basic  with  the  modified  vehicle  in  Tables  VI  and  VII 
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Table  VI 


Longitudinal  FLEXSTAB  Derivatives 


Derivative 


Hand  Calculated 
(rad”^) 


FLEXSTAB  Basic 
( rad“^ ) 


FLEXSTAB  SFG 
(rad~^) 


.2939 


.7448 


.1188 


.2152 


.0165 


.0143 


1.0772 


.1346 


.0736 


1.1158 


.1371 


.0740 


4.933 


5.0585 


5.4472 


-1.115 


-1.6566 


-1.7219 


6.4916 


10.4086 


10.6842 


.6559 


-17.1028 


-19.119 


-19.1128 


.0696 


.0762 


-2.513 


-2.696 


-2.682 


Table  VII 


Lateral  FLEXSTAB 

Derivatives 

Derivative 

Hand  Calculated 
(rad~^) 

FLEXSTAB  Basic 
(rad“^) 

FLEXSTAB  SFG 
(rad~^^) 

C 

-.4122 

-.4136 

-1.3428 

-.0418 

-.0434 

- .0723 

3 

C 

.068 

.0006 

.0158 

"3 

-.011 

-.1002 

- .1546 

•^p 

Jtp 

-.452 

-.5225 

- .6412 

c 

.0052 

.0227 

.0082 

n 

P 

C 

.3713 

.4105 

.5225 

^r 

c. 

.1608 

.1495 

.1784 

r 

C 

-.2177 

-.1899 

- .1968 

n 

r 

C 

.0111 

.0237 

y6A 

c. 

.2458 

-.2376* 

- .2824* 

^6  A 

C 

-.0252 

-.0167 

- .0155 

"6  A 

C 

.1971 

.2584 

.2600 

y6R 

.0017 

.0006 

.0003 

R 

C 

-.1203 

-.1355 

- .1355 

"x  r, 

6 R 


* Note:  FLEXSTAB  defines  its  control  derivatives  as  negatives. 
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Table  VIII 


Modified  FLEXSTAB  LATERAL  Derivatives* 


Deriv- 

ative 

Hand 
Calc 'd 

FLEXSTAB  BASIC 

__  Base-  , 

•®^yFUS  line  ^‘^^yVT 

% 

•®^yFus 

FLEXSTAB 

Base- 

line 

SFG 

C 

^3 

-.4122 

-.3784 

-.4136 

-.5875 

-1.2887 

-1.3428 

-1.5179 

-.0418 

-.0272 

-.0434 

-.0446 

- .0556 

- .0723 

- .0734 

3 

.068 

.0576 

.0006 

.0837 

.0707 

.0158 

.0987 

C 

-.011 

-.0994 

-.1002 

-.1122 

-.1532 

-.1546 

- .1665 

C 

s 

-.452 

-.5222 

-.5225 

-.5226 

-.640 

-.6412 

- .6414 

c 

n 

P 

.0052 

.0238 

.0227 

.0288 

.0092 

.0082 

.0142 

C 

^r 

.3713 

.4053 

.4105 

.6062 

.5141 

.5225 

.7194 

r 

.1608 

.1482 

.1495 

.1508 

.1771 

.1784 

.1795 

C 

n 

r 

-.2177 

-.1942 

-.1899 

-.2846 

-.2008 

-.1968 

-.2912 

^y6R 

.1971 

.2584 

.2584 

.3876 

.2600 

.2600 

.3900 

.0017 

.0006 

.0006 

.0009 

.0003 

.0003 

.0006 

c 

6R 

-.1203 

-.1355 

-.1355 

-.1891 

-.1353 

-.1353 

-.2031 

* Based  on  the  flight  condition  in  Table  II. 
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flight  test  the  vehicle  it  was  theorized  that  the  RPV  does  have  good 

weather  cock  stability.  Considering  this,  attempts  were  made  to  match 

the  FLEXSTAB  values  closer  to  the  hand  calculated  figures.  First  the 

influence  of  the  forward  vertical  portion  of  the  fuselage  was  reduced 

by  50%.  This  alteration  mathematically  states  that  there  is  less 

vertical  surface  area  forward  of  the  CG  which  increases  the  term. 

B 

The  second  modification,  increasing  the  effect  of  the  vertical  tail 
by  50%,  was  prompted  because  the  tail  section  experiences  a higher 
dynamic  pressure  due  to  the  prop  wash.  Neither  of  these  alterations 
had  any  effect  on  the  longitudinal  or  aileron  control  characteristics 
of  the  vehicle.  Table  VIII  displays  the  12  derivatives  that  were 
changed.  The  C terms  corresponding  to  the  altered  AIC  matrices 

"b 

(.5C  + 1.5C  ) agree  much  closer  to  the  hand  calculated  data. 

^FUS  ^VT 

Another  useful  set  of  data  obtained  from  the  SD  £ SS  program  is 
the  list  of  static  stability  parameters  in  Table  IX.  The  values  are 
given  in  decimal  fractions  of  the  mean  aerodynamic  chord,  c. 


Table  IX 

Static  Stability  Parameters* 


FLEXSTAB 

BASIC 

FLEXSTAB 

SFG 

Static  Margin 

.3161 

.3275 

Neutral  Point 

.6122 

.6175 

Load  Factor  (g) 

1.0022 

1.0021 

Maneuver  Point 

.6887 

.6997 

Maneuver  Margin 

.3926 

.4098 

Longitudinal  Control  per  G 

Steady  Pullup  (Deg/unit  load  factor) 

-4.9901 

-4.9953 

Stick  Speed  Stability  (Deg/in/sec) 

.0064 

.0064 

* Based  on  flight  condition  in  Table  I-I. 
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From  these  parameters  one  can  see  that  the  fore/aft  position  of  the 
CG  can  be  moved  aft  to  some  degree  without  making  the  vehicle  statically 
unstable  (CG  aft  of  Neutral  Point)  (Ref  10:209).  The  last  area  of  infor- 
mation gained  from  the  SD  S SS  program  is  the  characteristic  equation 
rooting  analysis.  As  the  name  infers,  the  characteristic  equation  of 
the  vehicle  is  analyzed  both  for  the  longitudinal  and  lateral  modes. 

Table  X displays  the  two  oscillatory  longitudinal  modes,  short  period 
and  phugoid. 


Table  X 

* 

Longitudinal  Modes* ** 


Model 

Short 

Period  (sec) 

Period 

Damping  Ratio 

r • • 

Phugoid 

Period  (sec)  Damping  Ratio 

1 

Basic 

1.28 

.5228 

15.2 

.0515 

SFG 

1.25 

.517 

15.2 

.0499 

* Note:  The  longitudinal  model  analysis  is  unchanged  by  the  alteration 
of  the  AIC  matrices. 


**  Base  of  flight  condition  in  Table  II. 

It  is  evident  from  Table  X that  only  a slight  variance  was  noted  in 
the  longitudinal  model  values.  The  SFG  modification  caused  a 2.3% 
decrease  in  the  short  period  period  and  a 1%  and  3%  decrease  in  the 
damping  ratios  of  the  short  period  and  phugoid  respectively. 

Table  XI  lists  the  lateral  characteristic  modes  which  are  the 
oscillatory  dutch  roll  mode  and  the  remaining  two  non-oscillatory  modes, 
rolling  mode  and  spiral  stability. 

It  is  evident  from  this  table  that  there  is  a significant  difference 
between  the  data  from  the  Baseline  AIG  matrix  and  the  two  altered 


Table  Xl 


Lateral  Modes* 


Roll 

Spiral 

Model 

J 

Dutch 

Roll 

Time  to 

Time  to  '4/2 

Period  (sec) 

Damping  Ratio 

Amplitude  (sec) 

Amplitude  (sec) 

Basic 

*^^FUS 

2.88 

.3834 

.078 

(2A) 

14.86* 

Baseline 

12.6 

.1578 

.078 

C4A) 

.427 

1.5C„ 

■'VT 

2.41 

.4144 

.078 

(2A) 

42.13** 

SFG 

2.63 

.4004 

.091 

(2A) 

29.29** 

•'FUS 

r4A) 

Baseline 

5.92 

.6680 

.091 

3.499 

l.SCv 

2.25 

.4238 

.091 

(2A) 

397** 

Yvt 

• 

- 

* Based  on  flight  condition  in  Table  II. 

**  Note:  These  spiral  modes  are  slightly  unstable. 


matrices.  Since  the  last  two  sets  of  data  were  obtained  by  using 
derivatives  that  were  closest  to  the  hand  calculated  values,  they  will 
be  discussed  further. 

The  dutch  roll  (DR)  period  decreased  approximately  8%  while  the  DR 
damping  increased  about  4%.  The  roll  response  was  decreased  because  the 
time  to  half  amplitude  for  the  roll  mode  was  increased  by  16.7%. 

The  last  lateral  mode  to  be  considered  was  the  spiral  mode.  The 
predominant  factor  in  determining  spiral  stability  is  the  magnitude 
and  sign  of  the  constant  term  in  the  characteristic  equation.  This 
constant  contains  the  following  combination  of  derivatives: 
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If  its  sign  is  positive  *he  vehicle  usvially  possesses  spiral  stability, 
but  if  the  value  is  negative  the  vehicle  will  display  spiral  instability. 
The  Baseline  model  has  spiral  stability  while  the  other  models  did  not. 
This  slight  spiral  instability  is  of  little  concern,  in  fact,  several 
general  aviation  aircraft  exhibit  an  unstable  spiral  mode  of  this 
magnitude  and  still  remain  easy  to  manually  control  See  Appendix  B 
for  a flow  chart  of  the  SD  S SS  program. 

Time  Histories 

The  Time  Histories  (TH)  program  calculates  the  nonlinear  time 
responses  for  coupled  perturbations  due  to  control  inputs,  discrete 
wind  gusts,  and/or  initial  perturbations.  The  TH  program  uses  infor- 
mation generated  by  the  GD  and  SD  S SS  programs  to  calculate  nonlinear 
coupled  perturbations  about  steady  state  conditions  in  the  time  domain. 
The  program  uses  this  inf ormation  to  generate  the  equations  of  motion 
which  include  gust  and  control  surface  forcing  function  coefficients. 

The  next  chapter  provides  the  detailed  time  histories  analysis.  See 
Appendix  B for  a flow  chart  of  the  TH  program. 
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perturbation  quantities  being  zero.  From  this  input  yaw  rate  (r), 
roll/yaw  coupling  and  side  slip  were  investigated  using  the  1.5C^ 
models. 

Figure  13  compares  the  yaw  rate  (r)  of  both  vehicles.  The  modified 
vehicle  reached  a yaw  rate  of  -.2061  rad/sec  at  five  seconds  compared  to 
the  basic  model's  yaw  rate  of  -.1864  rad/sec.  That  represents  a 10.6% 
increase  in  yaw  rate  due  to  the  SFG.  It  is  also  evident  that  the  basic 
vehicle  has  more  overshoot  (less  dutch  roll  damping). 

A comparison  of  the  roll  (Euler)  angle  phi  in  Figure  14  shows  that 
the  basic  model  rolled  through  -35.68*  while  the  SFG  model  achieved  a 
bank  angle  of  -37.03  , a 1.4%  increase  in  roll  angle  due  to  yaw. 


This  roll  due  to  yaw,  caused  by  the  increased  C.  term,  is  usually 

considered  an  adverse  effect  and  a possible  solution  to  the  problem 
would  be  to  lower  the  centers  of  pressure  of  the  SFG  by  increasing  their 
respective  spans.  Although  this  violates  an  earlier  restriction  from 
Chapter  II,  the  AFFDL  is  considering  placing  small  wheels  on  the 
extended  SFG  and  evaluating  this  option  in  the  wind  tunnel. 

Figure  15  is  a representation  of  the  side  slip  angles  beta.  The 
unmodified  vehicle  developed  a side  slip  angle  of  2.75*  in  five  seconds 
while  beta  of  the  SFG  model  was  2.06*,  a 25.1%  decrease.  The  SFG  vehicle 

possesses  more  yaw  stiffness,  mainly  due  to  the  increase  in  C , and  this 

"0 

accounts  for  the  smaller  side  slip  angle.  The  comp  time  history  response 
caused  by  the  rudder  step  is  located  in  Appendix  C. 


o 

Figure  15.  Side  Slip  Angle  Due  to  2 
Rudder  Step  (0  to  5 sec) 
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Aileron  Step 

An  aileron  step  of  two  degrees  was  applied  to  gain  knowledge  about 
the  vehicles'  roll  characteristics.  The  control  input  was  held  for  five 
seconds  starting  at  time  equal  zero  and  the  equations  were  integrated 
for  twenty  seconds.  The  reader  should  be  aware  of  the  sign  convention 
employed  by  FLEXSTAB;  all  the  control  derivatives  are  negative 
(Ct  , C , C < 0).  That  is,  a positive  aileron  deflection  will 

6a  "5r  “dE 

cause  a negative  roll.  Some  authors,  such  as  Seckel  (Ref  4)  and  Roskam 

(Ref  6),  do  not  follow  this  convention  and  define  C as  a positive 

’'6A 

quantity. 

Figure  16  is  a comparison  of  the  bank  angles,  phi,  for  the  1.5  C 

O ^VT 

vehicles.  At  the  termination  of  the  2 aileron  step,  the  basic  model 
achieved  a bank  angle  of  -73.3*  while  the  SFG  vehicle  rolled  through 

a -68.4*  angle. 

1.5  C„  Vehicles 
XVT 


The  7.2%  decrease  in  roll  ii-.3ponse  can  be  mainly  attributed  to  the 
increase  in  roll  damping,  C . This  loss  of  roll  response  was  expected 

because  when  any  lifting  (vertical  or  horizontal)  surfaces  are  added  to 
a vehicle  the  magnitude  of  the  roll  damping  increases. 

The  effect  of  this  roll"  damping  is  also  evident  in  the  roll  rate 
comparisons  of  Figure  17. 


l.^CyvT  Vehicles 


D.  r. 

11 

Basle  Vehicle 

20. 

• Solid  Line 

TIME 

SECONOS 

SFG  Vehicle 

• Broken  Line 

Figure  17.  Roll  Rate  Due  to  2* 
Aileron  Step  (0  to  5 sec) 


The  basic  vehicle  reaches  a value  of  -.290  rad/sec  and  the  SFG  model 
achieved  -.259  rad/sec  at  5 seconds.  This  is  a decrease  of  10.7%  in 
roll  rate  response. 

Another  factor  that  effects  roll  response  is  the  rolling  moment  of 
inertia,  I^x*  modification  increased  the  Ij^j^  value  by  43% 

which  means  theoretically  that  the  modified  vehicle  should  have  a lower 
roll  acceleration.  There  i.s  no  visible  change  in  the  slopes  of  the  roll 
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rate  curves  of  Figure  17  in  the  acceleration  periods  from  0 to  0.5  sec 
and  5.0  to  5.5  sec.  A conclusion  drawn  from  Figure  17  is  that  the 
increased  term  has  little  effect  on  the  roll  response. 

It  should  also  be  noted  that  both  airplanes  have  excellent  transient 
characteristics  in  that  they  have  little  overshoot  and  reach  their  steady 
state  roll  rate  (±  .01  rad/sec)  in  approximately  a half  second. 

The  last  area  considered  in  the  aileron  step  situation  was  adverse 
yaw  which  is  primarily  influenced  by  the  term.  Figure  18,  yaw  rate 

due  to  an  aileron  step,  shows  the  adverse  yaw  in  the  first  second  after 
control  application. 


Figure  18.  Yaw  Rate  Due  to  2* 
Aileron  Step  (0  to  5 sec) 


The  yaw  rates  generated  at  t - .3  sec  arc  .0074  rad/scc  and  .0068  rad/sec 


about  8.1%  more  adverse  yaw  rate  than  the  SFC  model,  but  considering  the 
very  small  magnitude  of  both  values,  adverse  yaw  should  not  be  a problem 
with  either  vehicle.  The  complete  time  history  response  caused  by  the 


aileron  step  is  located  in  Appendix  0. 

Cross  Control  Step  Input 

With  the  correct  combination  of  rudder  and  aileron  the  airplane 
will  execute  a flat  or  skidding  turn  with  a bank  angle  of  zero  degrees. 

To  determine  the  proper  control  deflections  to  accomplish  this  maneuver 
the  steady  state  problem  was  analyzed  by  representing  the  lateral 
equations  of  motion  in  state  variable  form: 

X = A X + B u (8) 

where  x is  a column  matrix  of  state  variables  B,  p,  r,  and  ♦.  u is  a 
column  matrix  containing  the  control  variables,  5 A and  5r.  A and  B are 
the  coefficient  matrices  containing  the  stability  and  control  derivatives 
respectively. 

Assuming  that  p,  and  the  rates  of  change  of  the  state  variables 
are  zero,  the  equations  can  be  written  in  the  following  form: 


0 = Cy 
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Using  a given  SR  of  .1  radians,  equations  (9),  (10),  and  (11)  were 

solved  simultaneously  for  ,the  remaining  3 unknown  variables.  From  this 

analysis  an  aileron  input  of  -.042  radians  would  re  required  to  maintain 

a flat  turn  of  zero  degrees  bank  with  a rudder  input  of  .1  radians  for 

the  basic  1.5  C model.  A ratio  of  this  combination  of  rudder  and 

y^T 

aileron  was  then  used  to  analyze  both  1.5  C vehicles. 

yyi 

Figures  19,  20,  21,  and  22  display  the  respective  comparisons  of  roll 
rate,  yaw  rate,  roll  angle  and  side  slip  angle. 
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- Broken  Line 

Figure  19.  Roll  Rate  due  to  4 * Rudder/ 
-1.68*  Aileron  Step  (0  to  5 sec) 
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Figure  22.  Side  Slip  Angle  Due  to  4® 
Rudder/-! .68®  Aileron  Step  (0  to  5 sec) 


From  these  figures  one  can  see  that  the  steady  state  control  input 

prediction  for  the  basic  1.5  C vehicle  was  in  error.  The  basic 

Yvt  ^ 

model  rolled  through  an  angle  of  -7.13  during  the  control  application. 
The  reason  the  steady  state  solution  differs  from  the  FLEXSTAB  solution 
is  because  FLEXSTAB  solves  the  entire  set  of  non-linear  differential 
equations  (Ref  6:2.24)  while  the  steady  state  method  uses  only  the 
uncoupled  linearized  lateral  equations  of  motion,  consequently  some 
small  but  significant  terms  are  neglected. 

One  method  of  overcoming  the  problem  of  proper-  control  inputs 
would  be  to  use  an  automatic  control  system  similar  to  the  one  depicted 
in  Figure  23.  With  this  controller  the  bank  angle  is  maintained  at  tlie 
specified  value  by  electronically  comparing  the  actual .bank  angle,  as 
sensed  by  the  gyro  system,  to  the  angle  desired.  The  autopilot  then 
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Figure  23.  Automatic  Control  System 


uses  this  terror  signal,  <|>  command,  to  generate  an  appropriate  aileron 
input  to  keep  the  error  signal  zero  and  the  bank  at  the  amount  specified. 

Notice  also  in  Figures  19  through  22  that  small  changes  in  the 
stability  derivatives  between  the  two  models  produced  a substantial 
change  in  the  rolling  and  yawing  characteristics.  A complete  time 
history  response  caused  by  the  cross  control  input  is  located  in 
Appendix  E. 


Gust  Inputs 

Two  wind  gusts  of  equal  magnitude  but  different  direction  were 
programmed  to  disturb  the  model  from  the  steady  state  flight  condition. 
In  the  first  case  a gust  of  11.09  ft/sec  from  the  positive  "y"  direction 
was  applied  to  the  model  for  one  second  starting  at  time  equal  zero. 
Refer  to  Figures  24,  25  and  26  for  roll  rate,  yaw  rate,  and  side  slip 
angle  histories,  respectively.  From  these  figures  it  is  easy  to  recog- 
nize the  oscillatory  dutch  roll  period  of  approximately  2.4  seconds  in 
the  two  to  five  second  intervals.  Figure  26  also  shows  that  the  SFG 
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1.5  Vehicles 


Figure  24.  Roll  Rate  Due  to  Side  Wind 
Gust,  11.09  fps  (0  to  1 sec) 
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Basic  Vchlcla  - Solid  Lina 
SFC  Vahlcla  - tirokan  Lina 


Figure  25.  Yaw  Rate  Due  to  Side  Wind 
Gust,  11.09  fps  (0  to  1 sec) 
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Figure  26.  Side  Slip  Angle  Due  to  Side 
Wind  Gust,  11.09  fps  (0  to  1 sec) 


vehicle  has  slightly  more  dutch  roll  damping  (amplitude  of  the  oscilla- 
tions are  lower)  and  a lower  natural  frequency. 

Because  the  six  equations  of  motion  are  coupled,  this  particular 
gust  also  excited  the  two  oscillatory  longitudinal  modes,  phugoid  and 
short  period.  The  phugoid,  characterized  by  a changing  forward  speed  u 
(Figure  27),  pitch  rate  q (Figure  28),  and  pitch  angle  0 (Figure  29), 
has  a period  of  approximately  fifteen  seconds  and  is  visible  after  five 
seconds  in  all  three  figures.  The  oscillatory  short  period,  as  seen  in 
the  1.5  to  4.0  second  time  interval  of  Figures  28,  29,  and  30,  is  well 
damped  and  has  a period  of  approximately  1.25  seconds. 
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Vehicles 


Basie  Vehicle  - Solid  Line 
SPG  Vehicle  - Broken  line 


Figure  28.  Pitch  Pate  Change  Due  to  Side 
Gust,  11.09  fps  (0  to  1 sec) 
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The  short  period  mode  is  again  observable  in  the  1.5  to  4.0  seeond 
intervals  of  Figures  32,  33,  and  35.  A complete  Lime  history  response 
due  to  the  gust  inputs  is  located  in  Appendix  F. 


Figure  34.  Angle  of  Attack  Due  to  Vertical 
Gust,  11.09  fps  (0  to  1 sec) 
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Basle  Vehicle.-  Solid  Lins 

XI  HE  SECONDS  Venlcls  - Broken  Lins  t 

Figure  35.  Vertical  Speed  Perturbation  I 

Due  to  Vertical  Gust,  11.09  fps  (0  to  1 sec)  j 
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VI.  Conclusions  and  Recommendations 

The  conclusions  drawn  in  this  section  were  based  on  the  time  history 
calculations  and  plots  generated  by  the  FLEXSTAB  program.  In  general  the 
modified  XBQM-106  displayed  an  overall  increase  in  lateral  stability  and 
response  while  suffering  a slight  reduction  of  longitudinal  stability. 

The  most  significant  effect  of  the  SFG  modification  was  to  increase 
the  yaw  rate  induced  by  a rudder  step  by  10.6%.  The  roll  due  to  yaw, 
an  adverse  condition,  was  increased  by  1.4%,  but  this  problem  could 
easily  be  overcome  by  a simple  autopilot. 

The  response  due  to  an  aileron  step  indicated  that  the  modification 
decreased  the  vehicle’s  roll  rate'  by  10.7%,  mainly  caused  by  the  increased 

term.  The  aileron  input  also  showed  the  vehicle's  adverse  yaw. 

P 

Although  these  yaw  rates  were  small  (a  maximum  of  less  than  0.4  deg/sec) 
the  basic  model  had  a 8.1%  higher  adverse  yaw  rate  when  compared  to  the 
SFG  vehicle. 

The  damping  ratios  of  both  longitudinal  modes,  short  period  and 
phugoid,  were  decreased  from  1%  to  3%  by  the  SFG  modification.  Examining 
the  lateral  modes,  the  dutch  roll  period  decreased  8%  while  the  DR  damping 
ratio  increased  by  4%.  The  roll  mode  time  to  half  amplitude  increased 
by  16.7%  and  the  spiral  stability  increased  for  all  the  models. 

The  reader  should  bear  in  mind  that  this  study  is  theoretical  in 
nature  and  a judicious  use  of  this  data  along  with  future  wind  tunnel  and 
flight  test  data  will  allow  the  AFFDL  to  make  a thorough  evaluition  of 
the  XBQM-106. 
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Appondlx  A 


Inortl.i/Contrr  o£  Gravity  Location  Pro/;ram 


Tr';T  ,Tfl'»El  = INPUT  , T Apr  ? = OUTPUT) 

^T-lEri5I0N  X ( I"!  ) fX  ( J 7c  ) ,7  (170)  f A (17C  ) , 3 ( 170)  , I'.T  (170) 

CALL  CFNAPAVd’C,  TOTUr  ,/> , H,  WT  , X , Y , Z , XCG , YCG , ZCC) 

«;ttp 

FUT 

SM^POUTIUF  ~FN"tPAV  (T>)“.  ,TOTVIT  I A ,P.,WT  ,X  ,Y  , 7 , XCG  , YCGf733) 

PFAL  ‘•10'1X,"0*IY,>''0MZ,  IA?r,IX>lTf  ITiIXZ 

oi'lEUSrON  X (IU‘I)  , Y ( IN!' ) ,Z  (INM)  I A ( lUM)  ,3  (IMM)  , WT  (INN) 

TOTMT  = YOMX  = 'tOMY=HOM7=0 . 0 
UPtTE(2,3007) 

9007  POPNATC*-  TTEN  WEIGHT  X T 

">0  qc08  I = 1,IM'1 

REAO  (1  .qoor,)  A(  I)  ,0(1)  ,»'T  (I),X(I)  ,Y(I),7(I) 

933?  rOJYAT  (2A!' ,^G10.3)  ' 

9003  WPTTE  (?,o30C)  A ( t ) , a ( j ) , KT ( I ) ,X { I) , Y (I) ,7 ( I ) 

93C0  (”  ",2110,^013.3) 

no  900  9 1 = 1,  I'JH  " . 

MOMXsMOMX  4-  WT(I)*X(T)  • 

Ma4Y  = HO‘1Y  Wf(T)»Y(r) 

Mn'47=H0'17  + WT{I)*7(I)  • . 

9309  nTWT  = TDTwr  f WT(T) 

WRTTE  (2,9007)  TOTWT 

9333  FO^-IAT  ("/.M°ORAPT  WEIGVT  IS",G10 . 3,"OUNCES.") 

XOO=HOMX/TOTMT 

YCr,  = MOMY/TOTWT  ■ . 

7r;';  = HOHZ/TOTWT  . 

WRITE  (2,90':<v)  XCG,YCG,'’CG 

900A  FORMAT ("/AIRCRAFT  CG  IS  LOCATED  AT  X=", GIO . 3, ", Y=",  313 . 3 , ", Z = ”, 
1G13.3,".") 

IX=IY=IZ=TX7=0.0 
PO  900  1 ! = 1,I'IM 
MA9S=(WT  (I)/l  j.  0/32.  17'! 
ixY=ixz--iAS<;*  (X  (i)-yci'.)  ■ (7(i)-7CG) 

IX=TX  ♦ lA'Jr,  ( ( (Y  ( I) -YCr  ) ■*  (Y  ( I) -YCG)  ) + { (Z(  T ) -7CG)  • ( 7 (I)  - ZCG) ) ) 

Tv=i  Y ♦ YAss''((  (X(r)-xcr)*(X(i)-xcG))  ♦•(  (Z(i)  -7cg)  * (Z  ( i ) -zee.) ) ) 

9001  I''=IZ  ♦ MAn'M  ( (X(  I) -xro  ♦(X  (I) -XCG)  ) +(  ( Y(  r ) -YCG)  •(  Y (I ) --YCC)  ) ) 

write  (2,P0i''’)  :x,IY,I7,ixZ 

9032  FORMAT  ("/RnMCUTS  OF  rcr'^TlA  IN  SLUG-INCHES  SPUARED  ARE  I-XX  = " 
IGl ) . 3,",  I-YY=",r.lO.  3,"I-7Z=",G10..7,"I-X7=",G10.3) 

FEIURN 

END 
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Appontnx  B 


FLKXSTAB  Flowcharts 


This  section  contains  the  basic  riowcharts  for  the  FLKXSTAB  programs 


used  in  this  study.  The  common  computer  symbols  are  defined  below. 

Flowchart  symbols: 

n Precompiler  statements 

a set  of  statements  producing 
U FORTRAN  source  code. 

I I Computational  processing 

a block  or  portion  of  a program  or  subproi'ram  that  processes 
a set  of  instructions  or  equations. 

Input/output  processing 

a block  or  portion  of  a program  or  subprogram  that  processes 
input/output  requests. 


Decision 

a point  in  the  program  where  a branch  to  a’: : mate  paths  is 
possible  depending  upon  variable  conditions. 

Terminal 

the  beginning  or  end  of  a program  or  subprogram. 


CZ3 

Q 

Q 


Flow  direction 

the  direction  of  processing  or  data  flow. 


Offpage  connector 


Punched  cards 


Output  listings 


Magnetic  tape 

used  for  permanent  storage  of  data  and  pro;'r.iiiis. 


Disk 

used  for  temporary  storage 
of  data  and  programs. 
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Bl,  Airplane  Definition  Suction 


Overlay  (2,0) 

Main  l*i(i)-iain  ODIM  OT 

Piiipow  To  iiioniior  .nul  eDniiol  eweiiiioii  llirdtii'-h  llio 
conlip.tii'alion  plnl  iiuxlule 


Subroutine  LOCATE 


Subroutine  CONTJ^L 
(See  flowchart  4.  3-10) 


Subroutine  LOCATE 


F’.rurc  B3.  Geometry  Definition  Program  Flowchart 


rnp 
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Flc^jrc  B5.  MG  Correction  Matrix  Program  Flowchart 


BBi  AVAIUBlt  COPY 


Oveilay  (0,0) 

Main  l*roi'r;uu  SPSS 

Purpose- -To  iniliato  oxoctilinti  in  ilio  SDASS  pioj'.rain 


Start 


Subroutine  IHITAL  Initialize  parameters  and 

(See  flowchart  12.3-2)  Inter-moduTe  labeled 

common  blocks. 


Subroutine  I0L'.\rS 


Subroutine  CLP.! A3 


Assign  numbers  to  the  I/O 
units  and  initialize  the 

matrix  file  position 
indicators • 


Clear  automatic  matrix 
management  package 
(AMMP)  tables 


Subroutine  VCONfL 


Print  verlson 
control  header 


Flf^urc  R6.  SU  fi  SS  Profjram  Flowchart 
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Figure  C2,  Side  Speed  Perturbation  Due  to  2 Ruddrr  Stepi 
(0  to  5 sec) 


C Vehicles 

'KUS 


0.  i. 

to. 

19. 

Male  Vch*cU 

to. 

- fT  ilM  Un« 

TIME 

SECONDS 

Sre  VeMO 

• r-r-jk^n  Line 

Figure  C3.  Townward  Speed  Perturbation  Due  to  2 liidderStep 
(0  to  5 sec) 
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Figure  G12.  Forward  Speea  Perturbation  Due  to  2 Rudder  Step 
(0  to  5 sec) 
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Figure  C13t  Side  Speed  Perturbation  Due  to  2°  Rudder  Step 
(O  to  5 sec) 
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Figure  C18,  Pitch  Angle  Due  to  2 Rudder  Step  (0  to  5 sec) 
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Figure  C19.  Roll  Angle  Due  to  2 Rudder  Step  (O  to  5 sec) 
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i*’!;;ure  C2l.  An/^le  of  Attack  Due  to  2 Hudder  Ste;;  (O  to  5 sec) 
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Figure  C24,  Side  Speed  Perturbation  Dus -to  2^  Rudder  Step 
(0  to  5 sec) 
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Figure  D4,  Roll  Rate  Due  to  2 Aileron  Step  (0  to  5 Sec) 
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Figure  D9.  Heading  Anf'le  Due  to  2®  Aileron  Step  (0  to  5 sec) 
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Figure  D14,  Downward  Speed  Perturbation  Due  to  2^Aileron  Step 
(0  to  5 sec) 


Figure  D15t  Poll  Rate  Due  to  2*  Aileron  Step  (O  to  5 -cc) 
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Figure  Dlfe.  Pitch  Rate  Due  to  2*  Aileron  Step  (O  to  5 sec) 


Figure  Dl?,  Yaw  Rate  Due  to  2**  Aileron  Step  (0  to  5 sec) 
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Figure  D24,  Side  Speed  Perturbition  Due  to  2 Aileron  Step 
(0  to  5 sec) 
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Fipjure  T)26,  Holl  Rate  Due  to  2*  Aileron  Step  (O  to  5 sec) 
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Appendix  E 

Time  History  Plots  for  Gross  Control  Step  Input 
This  section  contains  the  time  history  plots  for  the  modified  vehicles. 
The  Rudder/-!  .68 **  Aileron  step  Input  was  Initiated  at  time  zero  and 
held  for  5 seconds.  The  following  variables  are  plotted  In  this  order j 
u,  v,w,  p,Q 
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Figure  El,  Forward  Speed  Perturbation  Due  to  4^  Rudder/-! .68* 
Aileron  Step  (0  to  5 sec) 
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Figure  E4,  Roll  Rate  Due  to  4 Rudder/-!  .68*’  Aileron  Step 
(0  to  5 sec) 
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Figure  E6.  Yaw  Rate  Due  to  4' Rudder/-! .68*  Aileron  Step 
(O  to  5 sec) 


Figure  ElO.  Angle  of  Attack  Due  to  4 Rudder/-! ,68^  Aileron 
Ste  (0  to  5 sec) 
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Figure  Ell.  Side  Slip  Angle  Due  to  4*  Rudder/-!, ob  Aileron 
Step  (0  to  5 sec) 
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Figure  E12.  Forward  Speed  Perturbation  Due  to  4 Rudder/-!. 68 
Aileron  Step  (0  to  5 sec) 
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Figure  E13*  Side  Speed  Perturbation  Due  tC  4*  Rudder/-!  .68  ** 
Aileron  Stop  (O  to  5 sec) 
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Figure  El6.  Pitch  Rate  Due  to  4 Rudder/-! .68  Aileron  Step 
( 0 to  5 sec) 
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Figure  E17.  Yaw  Rate  Due -to  4^  Rudder/-! ,68*  Aileron  Step 
(0  to  5 sec) 


Appendix  F 

Gust  Input  Time  Response  Plots 

This  section  contains  the  complete  set  of  time  history  plots  for 
the  baseline  and  modified  vehicles.  Wind  gusts  from  the  u,v,h  directions 
with  a magnitude  of  11.09  fps  were  applied  at  time  zero  and  maintained 
for  one  second.  The  following  variables  are  plotted  in  this  orden 
u,  v,w,p,a,r,0,^, 


Figure  FI,  Forv.ard  Speed  Perturbation  Due  to  a Forward  Gust 
11.09  fps  (0  to  1 sec) 
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Figure  F2.  Downward  Speed  Perturbation  Due  to^a  Forward  Gust 
11,09  fps  (0  to  1 sec) 
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Figure  Fib.  Side  Slip  Angle  Due  to  a Side  Gust  11.09  fps 
(O  to  1 sec) 
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Figure  FI 8,  Downward  Speed  Perturbation  Due  to  a Downward 
Gust  11,09  fps  (0  to  1 sec) 
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Figure  F20,  Pitch  Angle  Due  to  a Downward  Gust  11.09  fps 
( 0 to  1 sec) 
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H i//urn  F21.  Angle  of  Attack  Due  to  a Downward  Gust  11.09fps 
(O  to  1 sec)  , 
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